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1.  INTRODUCTION 


The  propagation  of  Gaussian  beams  has  been  extensively  explored  in  the 
literature.  Most  of  the  early  literature  neglects  the  effects  of  beam 
truncation  by  an  aperture  with  the  result  that  a  beam  remains  Gaussian  as  it 
propagates. 1  ^  Often,  it  is  the  propagation  of  collimated  beams  that  is 
discussed, ^  ^  but  lately  the  emphasis  has  been  on  focused  beams.^  ^  With 
few  exceptions ^  ^  the  discussion  is  limited  to  aberration-free  circular 
beams . 


In  spite  of  the  vast  amount  of  literature  on  the  propagation  of  Gaussian 
beams,  there  is  little  in  terms  of  a  quantitative  comparison  with  the 
propagation  of  a  uniform  beam.  Even  when  results  for  both  types  of  beams  are 
occasionally  given, ^  8>10,H  tjiey  appear  in  a  normalized  form  in  such  a  way 
that  some  of  their  quantitative  difference  is  lost.  The  objective  of  this 
report  is  to  compare  quantitatively  the  effects  of  diffraction,  obscuration,  and 
aberrations  on  the  propagation  of  uniform  and  Gaussian  beams.  Equations  are 
derived  for  the  axial  and  transverse  irradiance  distributions  for  a  focused 
beam.  Equations  for  encircled-power  distributions  are  also  given.  For  a 
meaningful  comparison  between  a  uniform  and  a  Gaussian  beam,  the  total  power 
transmitted  by  the  aperture  is  kept  fixed,  regardless  of  the  value  of  the 
obscuration  or  the  nature  of  the  beam.  It  is  shown  that,  in  the  focal  plane, 
the  irradiance  at  the  focus  and  in  its  vicinity  is  smaller  for  a  Gaussian  beam 


^GaskiH,  J.  D.,  Linear  Systems,  Fourier  Transforms,  and  Optics,  John  Wiley 
and  Sons,  NY,  1978,  Section  10-7. 

^Siegman,  A.  E.,  An  Introduction  to  Lasers  and  Masers,  McGraw-Hill  Book 
Company,  NY,  1971,  Section  8-2. 

^Dickson,  L.  D.,  "Characteristics  of  a  Propagating  Gaussian  Beam," 

Appl.  Opt.,  9_,  1970,  pp.  1854-1861.  (This  paper  considers  the  effects  of 
aperture  truncation  and  derives  a  condition  under  which  they  may  be 
neglected. ) 

^Williams,  C.  S.,  "Gaussian  Beam  Formulas  from  Diffraction  Theory," 

Appl,  Opt.,  12 ,  1973,  pp.  871-876. 

^Herman,  R.  M. ,  J.  Pardo,  and  T.  A.  Wiggins,  "Diffraction  and  Focusing 
of  Gaussian  Beams,"  Appl .  Opt . ,  24 ,  1985,  pp.  1346-1354. 

^Buck,  A.  L. ,  "The  Radiation  Pattern  of  a  Truncated  Gaussian  Aperture 
Distribution,"  Pro c.  IEEE,  55 ,  1967,  pp.  448-450. 
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than  that  for  a  uniform  beam.  However,  away  from  the  focus  but  within  the 
central  disc  of  the  diffraction  pattern,  the  irradiance  is  higher  for  a  Gaussian 
beam.  Accordingly,  the  encircled  power  is  higher  for  a  uniform  beam  for  very 
small  circles,  but  the  reverse  is  true  for  larger  circles.  The  results  for  a 
collimated  beam  are  obtained  as  a  limiting  case  of  a  focused  beam,  namely,  a 
beam  focused  at  infinity.  The  differences  between  the  diffraction  patterns  for 
uniform  and  Gaussian  beams  decrease  as  the  obscuration  increases. 

The  problem  of  aberration  balancing  is  discussed  for  the  two  types  of 
beams.  Zernike  polynomials  representing  balanced  primary  aberration  for  uniform 
and  Gaussian  annular  beams  are  described.  The  relationship  between  the  peak 
value  of  a  primary  aberration  and  its  corresponding  standard  deviation  across  a 
circular  aperture  is  tabulated  for  both  uniform  and  Gaussian  beams.  It  is  shown 

that  a  Gaussian  beam  for  which  the  irradiance  at  the  edge  of  the  aperture  is 
_2 

e  times  the  irradiance  at  its  center  (if  there  were  no  obscuration)  is  only 

somewhat  less  sensitive  to  aberrations  than  a  corresponding  uniform  beam.  For  a 

Gaussian  beam  with  a  much  weaker  truncation  by  the  aperture,  the  aberration 

tolerance  increases  rapidly.  However,  in  that  case,  since  the  beam  power  is 

concentrated  in  a  small  region  near  the  center  of  the  aperture,  the  effect  of 

aberration  in  its  outer  region  is  negligible.  Accordingly,  for  a  weakly 

truncated  Gaussian  beam,  the  aberration  coefficients  are  defined  in  terms  of  the 
.  -2  . 

peak  aberration  at  the  e  irradiance  point  in  the  aperture  rather  than  at  its 
edge. 


^Campbell,  J.  P. ,  and  L.  G.  DeShazer,  "Near  Fields  of  Truncated  Gaussian 
Apertures,"  J.  Opt.  Soc.  Am.,  59 ,  1969,  pp.  1427-1429. 

®Olaofe,  G.  0.,  "Diffraction  by  Gaussian  Apertures,"  J.  Opt.  Soc.  Am.,  60, 
1970,  pp.  1654-1657. 

^Schell,  R.  G. ,  and  G.  Tyra,  "Irradiance  from  an  Aperture  with  Truncated- 
Gaussian  Field  Distribution,"  J.  Opt .  Soc .  Am. ,  61 ,  1971,  pp.  31-35. 

l^Nayyar,  V.  P. ,  and  N.  K.  Verma,  "Diffraction  by  Truncated-Gaussian  Annular 
Apertures,"  J.  Optics.  9,  (Paris),  1978,  pp.  307-310. 

^Holmes,  D.  A.,  J.  E.  Korka,  and  P.  V.  Avizonis,  "Parametric  Study  of 
Apertured  Focused  Gaussian  Beams,"  Appl ■  Opt. ,  11 ,  1972,  pp.  565-574. 

^Li,  Y.,  and  E.  Wolf,  "Focal  Shift  in  Focused  Truncated  Gaussian  Beams," 
Opt.  Comm.,  42 ,  1982,  pp.  151-156. 

^Tanaka,  K. ,  N.  Saga,  and  K.  Hauchi,"  Focusing  of  a  Gaussian  Beam  Through  a 
Finite  Aperture  Lens,"  Appl.  Opt.,  24,  1985,  pp.  1098-1101. 


Note  that  the  diffraction  effects  of  a  uniform  or  a  Gaussian  beam  are 
equivalent  to  those  of  a  uniform  or  Gaussian-apodized  pupil,  respectively,  in  an 
imaging  system.  Thus,  for  example,  an  irradiance  distribution  transverse  to  the 
direction  of  beam  propagation  represents  a  corresponding  point-spread  function 
of  an  imaging  system. 


^Lowenthal,  D.  D.,  "Marechal  Intensity  Criteria  Modified  for  Gaussian 
Beams,"  Appl.  Opt.,  13_,  1974,  pp.  2126-2133,  2774. 

I’Lowenthal,  D.  D.,  "Far-Field  Diffraction  Patterns  for  Gaussian  Beams  in 
the  Presence  of  Small  Spherical  Aberrations,"  J.  Opt.  Soc.  Am.,  65 ,  1975, 
pp.  853-855. 

l^Sklar,  E.,  "Effects  of  Small  Rotationally  Symmetrical  Aberrations  on  the 
Irradiance  Spread  Function  of  a  System  with  Gaussian  Apodization  Over  the 
Pupil,"  J.  Opt.  Soc.  Am..  £5,  1975,  pp.  1520-1521. 
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2.  DIFFRACTION  EQUATIONS 


3.  APERTURE  DISTRIBUTION 


Consider  an  annular  aperture  of  inner  and  outer  radii  of  ea  and  a, 
respectively,  where  0  e  <  1  is  the  linear  obscuration  of  the  aperture.  As 
indicated  in  Figure  1,  a  beam  of  outer  and  inner  radii  a  and  ea ,  respectively, 
is  focused  at  a  distance  R  from  the  aperture  plane.  The  beam  is  aberration- 
free  when  a  spherical  wavefront  of  the  radius  of  curvature  R  is  centered  at 
the  point  of  observation  and  passes  through  the  center  of  the  aperture.  For  a 
fixed  total  power  P^  transmitted  by  the  aperture  regardless  of  the  value  of 
e  or  the  nature  of  illumination  (uniform  or  Gaussian),  the  irradiance 
distribution  at  the  aperture  for  uniform  and  Gaussian  illumination  may  be 
written 


I  (p)  *  (P  /A) /( 1  -  e2)  (4) 

u  o 


and 


I^(p)  =  (Pq/A)  f(y;e)  exp(-2yp2) 


(5) 


respectively.  Here 


A  =  Tia2  (6) 

is  the  area  of  the  unobscured  circular  aperture 

f(y;e)  =  2y/(.e~2ytZ-  e~2y )  (7) 

and  p  is  in  units  of  a  so  that  E  <  p  <  1,  The  parameter  y  characterizes 
the  truncation  of  a  Gaussian  beam  by  the  aperture.  If  we  define  a  Gaussian 

beam  radius  cu  as  the  radial  distance  from  the  beam  center  at  which  the 

-2 

irradiance  is  equal  to  e  of  the  value  at  the  center  (if  there  were  no 
obscuration),  then 

y  *  (a/u>)2  (8) 
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Figure  1.  Diffraction  Geometry 
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1. e.,  /y  represents  the  ratio  of  the  aperture  radius  and  Gaussian  beam 

radius.  Note  that  as  y  -*•  0,  f(y;£)  (1  ~  e2)  *  and  Eq.  (5)  reduces 

to  Eq.  (4).  Hence,  the  results  for  a  Gaussian  beam  reduce  to  those  for  a 
corresponding  uniform  beam  in  the  limit  y  0. 

It  is  evident  from  Eq.  (4)  that  the  aperture  irradiance  (in  units  of 

P  /A)  for  a  uniform  beam  increases  with  e  as  (1  -  e2)  ^ .  For  a  Gaussian  beam, 
o 

the  irradiance  decreases  exponentially  from  a  maximum  value  of  f  (Y;e)exp(-2y£:2) 
at  the  inner  edge  of  the  aperture  to  a  minimum  value  of  f (y;e)exp(-2y)  at  the 
outer  edge.  For  example,  as  illustrated  in  Figure  2a,  when  y  =  1  and  e  =  0, 
the  maximum  and  minimum  values  are  2.31  and  0.31,  respectively.  Thus,  the 
mirrors  in  a  high-power  Gaussian  beam  train  are  not  only  illuminated  unevenly, 
but  they  must  also  withstand  considerably  higher  irradiance  levels  compared  to 
those  for  a  uniform  beam.  For  an  annular  beam,  the  peak  irradiance  is  even 
higher,  as  illustrated  in  Figure  2b  for  y  -  1  and  e  =  0.5,  since  the  same 
total  power  is  now  distributed  across  a  smaller  area.  The  ratio  of  the  peak 
values  of  the  aperture  irradiance  for  Gaussian  and  uniform  beams  is  given  by 

I  (e)/I  ( e)  =  (1  -  e2)  f(Y;e)exp(-2ye2) 
g  u  ' 

=  2y(  1  -  G2 ) / { 1  ~  exp[-2Y(l  -  e2)]}  (9) 

The  variation  of  this  ratio  with  e  is  shown  in  Figure  2c  for  / y  =  0.5,  1, 

2,  and  3.  It  is  evident  that  as  yd  “  E2)  increases,  the  ratio  approaches 
a  value  of  2y(l  -  £2). 


-15- 


Several  Values  of  y 


4.  FOCUSED  BEAM 


For  an  aberration-free  beam  focused  on  the  axis  of  the  aperture  at  a 
distance  R  from  it,  the  amplitude  at  the  aperture  can  be  written 

U(p)  *  /iTpT  exp(-iAp^/XR)  (10) 

If,  for  example,  R,  is  the  radius  of  curvature  of  a  spherical  wave  incident  on 
b  -1-1-1 

a  lens  of  focal  length  f,  then  R  =  +  f  Substituting  Eq.  (10)  into 

Eq.  (1),  we  find  that  the  irradiance  distribution  in  an  observation  plane 
can  be  written 


I(r;z)  =  4(R/z) 


exp[  i<J>2(p)  ] 


J0(nrp) 


pdp 


(11) 


where 

-  x  (;  -  i}>2  112 

represents  the  defocus  phase  aberration  of  a  beam  focused  at  a  distance  R  with 
respect  to  a  reference  sphere  centered  at  a  distance  z  and  passing  through  the 
center  of  the  aperture.  In  Eq.  (11),  various  quantities  are  normalized  by  the 

parameters  for  a  uniformly  illuminated  circular  aperture.  The  I(p)  is  in  units 

2  2...  . 
of  P  /A,  I(r;z)  is  in  units  of  P  A/X  R  ,  p  is  in  units  of  a,  and  r  =  I r I 
o  o 

is  in  units  of  Xz/D,  where  D  *  2a. 


We  shall  refer  to  the  observation  plane  at  z  *  R  where  the  beam  is  focused 
at  as  the  focal  plane  and  the  axial  point  at  z  =  R  as  the  focus  or  the  focal 
point.  Similarly,  we  shall  refer  to  a  beam  with  R  =  00  at  the  exit  of  the 
aperture  as  a  collimated  beam.  In  imaging  applications,  the  focal  plane  will  be 
more  appropriately  called  an  image  plane,  unless  the  object  is  at  infinity. 
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4.1  FOCAL-PLANE  DISTRIBUTION 


If  we  let  z  =  R  and  substitute  aperture  distributions,  Eqs.  (4)  and  (5), 
for  uniform  and  Gaussian  beams  into  Eq.  (11),  we  obtain  the  corresponding 
focal-plane  distributions.  In  the  case  of  a  uniform  beam,  it  is  given  by^®*^ 


I  (r;R)  =  (1  -  e2)-1  [/TTrT  -  e2  /TTerT]2  (13) 

u  c  c 

where 

I  (r)  =  [2J1(itr)/irr]2  (14) 

C  l 

represents  the  focal-plane  irradiance  distribution  for  a  uniform  circular 
(e  =  0)  beam.  It  is  evident  that  for  the  focal-plane  distribution,  the 
variable  r  is  in  units  of  \F,  where  F  =  R/D  is  the  f-number  of  the  focused 
beam.  The  encircled  power,  in  units  of  PqI  is  given  by 

l 

P  (r  ;R)  =  (1  -  e2)  l[P  (r  )  +  e2P  (er  )  -  4  e  /  J,  (itrr  JJ^irerr  )dr/r]  (15) 

uo  co  co  ^lolo 

where  r  is  in  units  of  \F,  and 
o  ’ 

P  (r  )  -  1  -  J2(flr  )  -  J2(ttt  )  (16) 

co  0  o  l  o 

represents  the  focal-plane  enc ire led-power  distribution  in  the  case  of  a  uniform 
circular  beam. 


The  corresponding  results  for  a  Gaussian  beam  are  given  by 


V'!R) 


4f(y;e) 


JQ(T!rp)pdp 


2 


(17) 


l®Born.  M. ,  and  E.  Wolf,  Principles  of  Optics,  Pergamon  Press,  NY,  1975, 

p .  4 16 . 

i^Mahajan,  v.  N. ,  "Included  Power  for  Obscured  Circular  Pupils."  Appl.  Opt., 
_1_7 ,  197$,  pp.  964-968. 
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and 


r 


P  (r  ;R)  =  (it 
g  o’ 


o 

I  (r;R)rdr 
g 


(18) 


Figure  3  shows  how  the  irradiance  and  encircled-power  distributions  in  the 

focal  plane  compare  for  the  two  beams.  The  uniform  beam  corresponds  to  y  =  0, 

and  the  Gaussian  beam  considered  in  Figure  3  is  for  y  =  1.  At  and  near  the 

focal  point,  a  uniform  beam  gives  a  higher  irradiance  than  a  Gaussian  beam. 

For  a  circular  aperture  (Fig.  3a),  I  >  I  for  r  <  0.42.  For  larger  values 

u  g 

of  r,  I  >  I  ,  except  in  the  secondary  rings  where  again  I  >  I  ,  as  is 
gu’r  u  g 

well  known.  The  encircled  power  P  ^  P  for  r  $  0.63.  Of  course,  as  r 

u  v  g  o  o 

P  -*■  P  -*■  1.  The  Gaussian  illumination  broadens  the  central  disc  but  reduces 
u  g 

the  power  in  the  secondary  rings.  For  annular  apertures  with  e  =  0.5  (Fig.  3b), 

the  differences  between  the  focal-plane  distributions  for  uniform  and  Gaussian 

beams  are  less  compared  to  those  for  a  circular  aperture.  The  obscuration  reduces 

the  focal-point  irradiance,  reduces  the  power  in  the  central  disc,  and  spreads  it 

into  the  secondary  rings  of  the  diffraction  pattern.  It  reduces  the  size  of  the 

central  disc  also.  Moreover,  the  difference  in  encircled  powers  P  -  P 

u  g 

changes  its  sign  from  positive  to  negative  to  positive  as  rQ  increases.  Note 
also  that  because  of  the  obscuration,  the  secondary  maxima  are  higher  and  of 
nearly  equal  value  for  the  two  types  of  beams. 


For  clarity,  the  irradiance  distributions  are  also  plotted  on  a  logarithmic 
scale  (Figs.  3c  and  3d).  The  positions  of  maxima  and  minima  and  the  correspond¬ 
ing  irradiance  and  encircled-power  values  are  given  in  Table  1.  It  is  evident 
that  the  corresponding  maxima  and  minima  for  a  Gaussian  beam  are  located  at  higher 
values  of  r  than  those  for  a  uniform  beam.  Moreover,  whereas  the  principal 
maximum  for  a  Gaussian  beam  is  only  slightly  lower  (0.924  compared  to  1),  its 
secondary  maxima  are  lower  by  a  factor  >3  compared  to  the  corresponding  maxima 
for  a  uniform  beam. 


The  focal-point  irradiance  corresponding  to  uniform  and  Gaussian  beams  is 
given  by 


Iu(0;R)  =  1  -  E2 


(19) 
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CO  to  "T  C\i  o 
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Ere  ire led-Power  Distributions 
1)  Beams  of  a  Given  Total  Power 


21 


Table  1.  Maxima  and  Minima  of  Focal-Plane  Irradiance  Distribution  and 
Corresponding  Encircled  Powers  for  Uniform  and  Gaussian 
( Y  =  1)  Circular  Beams 


Maxima /Minima 

r>«'o 

I  ( r) 

P(r0) 

Max 

0 

1 

0 

(0) 

(0.924) 

(0) 

Min 

1.22 

0 

0.838 

(1.43) 

0 

(0.955) 

Max 

1.64 

0.0175 

0.867 

(1.79) 

(0.0044) 

(0.962) 

Min 

2.23 

0 

0.910 

(2.33) 

(0) 

(0.973) 

Max 

2.68 

0.0042 

0.922 

(2.76) 

(0.0012) 

(0.976) 

M  in 

3.24 

0 

0.938 

(3.30) 

(0) 

(0.981) 

Max 

3.70 

0.0016 

0.944 

(3.76) 

(0.0005) 

(0.983) 

Min 

4.24 

0 

0.952 

('  -29) 

(0) 

(0.985) 

Max 

4.71 

0.0008 

0.957 

(4.75) 

0.0002 

(0.986) 

Note:  The  numbers  without  parentheses  are  for  a  uniform  beam  and 
those  with  parentheses  are  for  a  Gaussian  (y  =  1)  beam. 
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and 


Ig(0;R)  =  (2/y)  t anh [ ( 1  -  e2)y/2]  (20) 

•  ,  .  .  .  .  u  20 
respectively.  Their  ratio  is  given  by 

0  =  I  (0;R)/I  (0;R) 
g  u 

=  tanhKl  -  e2)y/2)  (21^ 

(1  -  €2 )y/2 

Figure  4a  shows  how  r|  varies  with  y  for  several  values  of  e.  It  is 
evident  that  r|  decreases  as  y  increases,  regardless  of  the  value  of  e. 

However,  as  shown  in  Figure  4b,  for  a  given  value  of  y,  0  increases  as  e 
increases.  Note  that  for  large  values  of  y 

H  -*•  2/y ( 1  -  ez)  (22) 


The  decrease  in  r|  due  to  an  increase  in  y  can  be  due  to  an  increase  in 

a  and/or  a  decrease  in  uj.  If  we  consider  the  absolute  (unnormalized)  values  of 

the  focal-point  irradiance,  we  note  that  Iu(0;R)  increases  quadratically  with  a, 

but  Ig(0;R)  increases  only  slightly  (<_  10%)  with  a  for  /y  3.  Similarly, 

Ig(0;R)  decreases  (nearly)  quadratically  as  to  decreases  for  /y  _>  3.  It  is 

evident  that  n  is  always  less  than  1.  A  redistribution  of  the  aperture  power  from 

20 

a  uniform  to  any  nonuniform  distribution  reduces  the  focal-point  irradiance. 


4.2  AXIAL  IRRADIANCE 


If  we  let  r  =  0  in  Eq.  (11),  we  obtain  the  axial  irradiance  for  uniform 
and  Gaussian^  beams 


17 


I  (0;z)  *  (R/z)2(l  -  e2)  {sinf ( 1  -  e2)4>  /2]/[(l  -  e2>$  /2]}2  (23) 

u  o  o 


^^Maha  jan,  V.  N.,  "Luneburg  Apodization  Problem  I,"  Opt.  Lett.,  5,  1980, 
pp.  267-269. 
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and 


I  <0;z)  =  (R/z)2[2y/(<J>o2  +  y2)]{coth[(l  -  e2)y] 

-  cos  I  ( 1  -  e2)<t>  ]  /  sinh  [  ( 1  -  e2)y]}  (24) 

o 

respectively-  The  quantity  $0  represents  the  defocus  phase  aberration  at 
the  edge  of  a  circular  aperture.  It  is  given  by 


o 


<K  - 1> 


(25a) 


=  ttN(-  -  1) 
z 

where 


(25b) 


N  =  a2/XR 

is  the  Fresnel  number  of  the  circular  aperture  as  observed  from  the  focus. 


(26) 


The  positions  of  maxima  and  minima  of  axial  irradiance  are  obtained  by 
equating  to  zero  its  derivative  with  respect  to  z.  In  the  case  of  a  uniform 
beam,  the  minima  have  a  value  of  zero.  They  are  located  at  z  values 
corresponding  to  an  integral  number  of  waves  of  defocus  as  an  aberration  at 
the  outer  edge  of  the  aperture  relative  to  that  at  its  inner  edge,  i.e.,  for 


<J0  =  2im/(l  -  e2),  n  =  +1,  +2,  ...  (27a) 

or 

z/R  *  {1  ♦  [2n/N(l  -  e2)]}'1  (27b) 

The  positions  of  maxima  of  axial  irradiance  are  given  by  the  solutions  of 

tan((l  -  e2)4>0/2l  -  (R/z)(l  -  e2)4>0/2,  z  *  R  (28) 


In  the  case  of  a  Gaussian  beam,  the  positions  of  minima  and  maxima  of  axial 
irradiance  are  given  by  the  solutions  of 


{cosh ( ( 1 


e2)y] 


cos [ (l  -  e2 )$  ]} 
o 


*  -(1  -  e2)sin[(l  -  e2)$  ) 

o 


(24) 
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Figure  5  s'.ows  how  the  axial  irradiance  of  a  uniform  focused  beam  differs 
from  that  for  a  Gaussian  beam  when  y  -  1  and  Fresnel  number  N  =  1,  10,  and 
100.  We  note  that  the  principal  maximum  is  higher  for  the  uniform  beam 
compared  to  that  for  the  Gaussian  beam.  However,  the  secondary  maxima  are 
higher  for  the  Gaussian  beam.  Moreover,  whereas  the  axial  minima  for  the 
uniform  beam  have  a  value  of  zero,  the  minima  for  the  Gaussian  beam  have 
nonzero  values.  For  a  given  value  of  E,  the  locations  of  maxima  and  minima, 
except  the  principal  maximum,  are  very  nearly  the  same  for  the  two  beams.  The 
effect  of  the  obscuration  is  to  reduce  the  irradiance  at  the  principal  maximum 
but  increase  it  at  the  secondary  maxima.  Also,  the  maxima  and  minima  occur  at 
smaller  z  values  for  an  annular  aperture.  These  z  values  correspond  approxi¬ 
mately  to  those  axial  points  at  which  the  annular  aperture  subtends  an  odd  or 
an  even  number  of  Fresnel's  halfwave  zones,  respectively.  We  note  that  the 
curves  become  symmetric  about  the  focal  point  z  =  R  as  N  increases. 

Note  that  even  though  the  principal  maximum  of  axial  irradiance  does  not 

lie  at  the  focus,  maximum  central  irradiance  on  a  target  at  a  given  distance 

from  the  aperture  is  obtained  when  the  beam  is  focused  on  it.  This  can  be  seen 

by  equating  to  zero  the  derivative  of  axial  irradiance,  Eqs.  (23)  and  (24), 

with  respect  to  R.  When  doing  so,  the  normalization  factor  P  A/A2R2  should  be 

o 

substituted  in  these  equations  with  the  consequence  that  the  R2  factor  in 
front  of  their  right-hand  side  disappears.  Figure  6  illustrates  how  the 
central  irradiance  on  a  target  at  a  fixed  distance  z  varies  when  the  beam  is 
focused  at  various  distances  R  along  the  axis.  The  irradiance  in  this  figure 
is  in  units  of  PqA/A2z2.  The  quantity  N^  =  a  /Az  represents  the  Fresnel  number 
of  a  circular  aperture  as  observed  from  the  target.  As  in  Figure  2,  the  maximum 
irradiance  values  for  uniform  and  Gaussian  (y  =  1)  beams  are  1  and  0.924, 
respectively,  when  £  =  0,  and  0.750  and  0.717  when  £  =  0.5.  We  note  that  as 
increases,  the  curves  become  symmetric  about  R  =  z. 

It  is  evident  from  Eqs.  (23)  and  (24)  that  the  axial  irradiance  depends 
on  z  and  R  through  the  inverse-square-law  dependence  in  (R/z)  factor  and  the 
defocus  aberration  <5*^.  The  irradiance  is  symmetric  with  respect  to  the  sign 
of  <t>Q.  We  note  from  Figures  5  and  6  that  as  the  Fresnel  number  N  or  Nz 
becomes  large,  the  axial  irradiance  becomes  symmetric  about  the  point  z  =  R. 

The  reason  for  this  is  simple,  as  may  be  seen  by  an  examination  of  Eq.  (25). 
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Figure  5.  Axial  Irradiance  of  a  Beam  Focused  at  a  Fixed  Distance  R  with 
Fresnel  Number  N  =  1,  10,  and  100 
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Fipure  5 .  Axial  Trradiance  of  a  Ream  Focused  at  a  Fixed  Distance  R  with 
Fresnel  Number  N  =  1,  10,  and  100  (Continued) 


Figure  6.  Central  Irradiance  on  a  Target  at  Fixed  Distance  z  from  the  Aperture 
Plane  When  Beam  is  Focused  at  Various  Distances  R 


When  N  is  very  small,  z  must  be  much  different  from  R  for  4>  to  have  a 

o 

significant  value.  For  example,  when  N  m  1,  a  quarter  wave  of  defocus  = 

+  tt/2)  is  obtained  for  R/z  ■  1  ±  0.5.  Similarly,  when  N  *  1,  *  ±  tt/2  for 

z  o 

z/R  *  1  +  0.5.  In  either  case,  the  axial  irradiance  is  asymmetric  about  the 
point  z  -  R  because  of  the  asymmetry  in  the  defocus  aberration  as  well  as  in 
the  inverse-square-law  dependence.  However,  when  N  or  Nz  is  very  large,  a 
slight  difference  between  z  and  R  gives  a  large  value  of  For  example,  when 

N  ■  100,  a  quarter  wave  of  defocus  is  obtained  for  R/z  =  1  ±  0.005,  or  z/R  -  1 
+  0.005.  Note  that,  for  large  N,  |4>o  I  is  symmetric  about  the  focal  point. 
Hence,  in  the  vicinity  of  focus,  where  the  axial  irradiance  is  appreciable 
compared  to  the  focal-point  irradiance,  the  axial  irradiance  is  symmetric 
about  the  focal  point.  The  inverse-square-law  dependence  has  negligible 
effect  in  this  case,  since  z  and  R  are  practically  equal  to  each  other. 

Similar  considerations  hold  when  is  very  large,  i.e.,  the  irradiance  at 

the  point  of  observation  is  the  same  for  two  beams  focused  symmetrically  about 
it.  These  considerations  also  extend  to  the  transverse  irradiance  distribu¬ 
tions  (see  Eq.  (30)].  Thus,  for  example,  the  transverse  irradiance  distribu¬ 
tion  in  the  vicinity  of  the  focal  plane  is  symmetric  about  it  only  when  N  is 
very  large. 

4.3  DEFOCUSED  DISTRIBUTION 


When  z  t  R,  we  may  write  the  right-hand  side  of  Eq.  (11)  as  the 
product  of  two  integrals  and  retain  only  its  real  part,  since  irradiance  is  a 
real  quantity.  Thus,  the  irradiance  distribution  in  a  defocused  plane  can  be 
written 


I(r;z)  *  4(R/z)‘ 


// 


ps/TTpTiTsT  cos[4>  (p2 
o 


s2)  ]  Jgiirrp)  JQ(irrs)dpds 


(30) 


If  we  let  r  *  0  and  note  that  Jq(0)  =  1 ,  we  obtain  a  different  form  of  the 
expression  for  axial  irradiance,  namely 


l  .1 


I  (0;z)  =  4(R/z )‘ 


// 


ps/i(p)Ks) 


cos(<t  (p2  - 
o 


s2 ) )dpds 


(31) 
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The  encircled  power  (in  units  of  PQ  with  r0  in  units  of  Xz/D)  is  given  by 


P(r  ;z) 
o 

where ^ 

Q(p,s;r  ) 
o 

=  (r^/2)  [ J^(7Tr  p)  +  J^(Tir  s)]  if  p  *  s  (33a) 

o  0  o  lo 

*  [r  /ir(p2  -  s2)][pJ,(nr  p)  J-(irr  s) 

O  I  o  u  o 

-sJ,(rrr  s)  J-(7Tr  p)]  if  p  ?  s  (33b) 

1  o  0  o 


1  .1 


2tt "  J  J  ps/l(p)I(s)  cos[4>  (p2  -  s2)]Q(p,s;r  )dpds  (32) 
e  e  0  ° 


=  /  J0(7Trp)  J0 


(firs )rdr 


The  integrals  in  Eqs.  ( 30 ) — ( 32 )  can  be  evaluated  by  the  Gauss  quadrature  method 

7  7 

according  to  which  we  may  write"6 


1 


e  e 


f(p,s)  dp  ds  =  [(1  -  e)/2] ' 


y  u)Z  f(p .,s .) 
i-i  1  1  J 


M  i-1 

+  2  y  y  uj.u) .  f (p .  ,s . ) 

i=2  j=l  1  J  1  J 


(34) 


where  M  is  the  number  of  1-D  quadrature  points,  u).  are  the  weight  factors, 
and 


=  [1  +  e  +  (1  -  e)x{]/2  (35) 

x.  being  the  ith  zero  of  the  Mth-order  Legendre  polynomial.  In  our  calcula¬ 
tions  we  have  used  a  24-point  (M  *  24)  Gauss  quadrature.  Note  that  by  letting 

2 lAbramowitz ,  M.,  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions, 

Dover,  NY,  1970,  p.  484. 

22  Ibid  . ,  p.  887. 
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4>o  =  0  in  Eqs.  (30)-(32),  we  can  calculate  the  focal-plane  distributions  as 

well.  Equations  ( 30 ) — ( 35 )  are  generalizations  of  our  earlier  work  on  uniform 
.  23 

circular  beams  to  annular  apertures  and  Gaussian  beam.  Note  that  with 
slight  modification,  Eqs.  (30)-(32)  can  be  applied  to  diffraction  calculations 
involving  any  radially  symmetric  amplitude  and  phase  distributions  at  the  aper¬ 
ture.  For  example,  if  spherical  aberration  (A  p1*)  were  present,  the  cosine 

s 

factor  in  these  equations  would  become  cos[$(p)  -  $(s)]  where  $(p)  *  <t  p2  +  A  p1*. 

o  s 

An  example  of  defocused  distribution  is  illustrated  in  Figure  7  for  both 

uniform  and  Gaussian  (y  =  1)  beams  with  a  large  Fresnel  number  (so  that  the 

inverse-square-law  variation  is  negligible)  when  e  =  0  and  e  =  0.5.  The 

amount  of  defocus  $  =  2.783  rad  (or  0.443X)  is  such  that  the  central 

o 

irradiance  for  a  uniform  circular  beam  is  reduced  to  half  the  corresponding 

focal-point  irradiance.  (The  defocused  distributions  shown  can  also  be 

interpreted  as  the  distributions  on  a  target  at  a  fixed  distance  z  when  the 

beam  is  focused  at  a  distance  R  such  that  $  =  2.783  rad.  In  this  case, 

o 

the  irradiance  would  be  in  units  of  P  A/X2z2  and  r  and  r  would  be  in 

o  o 

units  of  Xz/D.)  We  note  that,  as  in  the  case  of  focal-plane  distributions, 
the  central  irradiance  for  a  Gaussian  beam  is  lower  than  that  for  a  correspond¬ 
ing  uniform  beam.  Note,  however,  that  the  defocus  aberration  does  not  reduce 
the  central  irradiance  for  the  annular  beam  as  much  as  it  does  for  the  circular 
beam,  so  much  so  that,  for  the  amount  of  defocus  aberration  considered  in 
Figure  7,  the  defocused  central  irradiance  for  the  annular  beam  is  higher  than 
that  for  the  corresponding  circular  beam.  For  the  uniform  and  Gaussian  circu¬ 
lar  beams,  the  central  irradiance  decreases  from  1  and  0.924  to  0.500  and  0.483, 
respectively.  For  the  annular  beams,  it  decreases  from  0.750  and  0.717  to 
0.514  and  0.497,  respectively.  This  indicates  the  well-known  fact  that  the 
tolerance  for  a  radially  symmetric  aberration  such  as  defocus  is  higher  for  an 
annular  beam  than  that  for  a  circular  beam.  When  E  =  0,  the  encircled  power 

is  higher  for  a  uniform  beam  for  small  values  of  r  compared  to  that  for  a 

o 

Gaussian  beam.  When  £  *  0.5,  the  difference  in  encircled  power  for  the  two 
types  of  beam  changes  from  positive  to  negative  to  positive  depending  on  the 
value  of  r  . 


^Mahajan,  V.  N. ,  "Aberrated  Point-Spread  Function  for  Rotationally  Symmetric 
Aberrations,"  Appl.  Opt.,  22,  1983,  pp.  3025-3041. 
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Figure  7.  Defocused  Irradiance  and  Encircled-Power  Distributions 
for  Uniform  and  Gaussian  (y  =  1)  Beams 


5.  COLLIMATED  BEAM 


The  results  for  a  collimated  beam  can  be  obtained  from  those  for  a  focused 
beam  by  letting  R  -*■  °°.  Thus,  for  example,  Eq.  (11)  for  the  irradiance  distri¬ 
bution  in  a  plane  at  a  distance  z  from  the  aperture  reduces  to 


where 


I  ( r ;  z ) 


/l(p)  exp[  i<t>o  p2  ] 


J0(iTrp)pdp 


(36) 


<t  =  A/Xz  (37) 

o 

represents  the  (negative  of)  phase  aberration  of  a  plane  wavefront  with 

respect  to  a  reference  sphere  centered  at  a  distance  z  and  passing  through  the 

center  of  the  aperture.  In  Eq .  (36),  the  irradiance  in  both  the  aperture  and 

the  observation  planes  is  in  units  of  the  aperture  irradiance  P  /A  for  a 

o 

uniform  circular  beam.  As  in  Eq .  (11),  r  is  in  units  of  Xz/D. 

In  the  far  field,  i.e.,  for  z  >  D2 /X ,  the  phase  aberration  $  <  tt/4 

—  o  — 

(corresponding  to  a  wave  aberration  of  less  than  or  equal  to  X/8)  and  may  be 
neglected.  Hence,  the  irradiance  distribution,  and  correspondingl  the 
enc ire led-power  distribution  in  a  far-field  plane,  is  similar  to  a  focal-plane 
distribution  discussed  earlier.  The  only  difference  is  in  scaling  of  the 
diffraction  pattern.  Similarly,  in  the  near-field,  i.e.,  for  z  <  D2/X,  the 
irradiance  and  encircled-power  distributions  correspond  to  defocused  distri¬ 
butions  discussed  earlier.  The  only  significant  difference  is  in  the  definition 

of  <t>  . 
o 

If  z  is  in  units  of  the  far-field  distance  D2 /X ,  and  we  let  r  =  0  in 

Eq.  (36),  we  obtain  the  axial  irradiance  (in  ,jnits  of  P  /A)  for  uniform17 

o 

and  Gaussian  beams 

lu(0;z)  -  (4/(1  -  e2)J  sin2[Tr(l  -  e2)/8z]  (38) 
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and 


l^(0;z)  =  {2y/[l  +  (Ayz/ir)2]}  {cothtyd  -  e2)] 

-  cos[tt(1  -  e2 ) /4z ] /  sinh  [y(  1  -  e2)]}  (39) 

respectively.  For  a  uniform  beam,  the  maxima  of  axial  irradiance  have  a  value 
of  4/(1  -  e2)  at 


z  =  (1  -  e2) /4( 2n  +  1)  ,  n  =  0,  1,  2,  ...  (40) 

Its  minima  have  a  value  of  zero  at 

z  =  (1  -  e2)/8n  ,  n  =  1,  2,  ...  (41) 

The  positions  of  maxima  and  minima  of  axial  irradiance  for  a  Gaussian  beam  are 
given  by  an  appropriate  modification  of  Eq .  (29),  namely 

{2(4z/tt)\2/[  1  +  (4yz/rr)2  ]  }  {cosh  [y(  1  -  e2)]  -  cos[it(1  -  e2)/4z)} 

*  -  (1  -  G2)sin[n(l  -  £2)/4z]  (42) 

Figure  8  illustrates  how  the  axial  irradiance  of  collimated  uniform  and 
Gaussian  beams  varies  with  distance  z  from  the  aperture.  With  reference  to 
Figure  5,  Figure  8  corresponds  to  N  -  0.  In  contrast  to  Figure  5,  the  maxima  of 
axial  irradiance  of  a  collimated  uniform  beam  have  the  same  value  of  4/(1  -  e2). 
Moreover,  unlike  the  principal  maximum  in  Figures  5a  and  5d ,  the  maximum 
farthest  from  the  aperture  has  a  lower  value  than  the  maxima  closer  to  it  in  the 
case  of  a  Gaussian  beam. 
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6.  ABERRATION  BALANCING 


A  focused  beam  emerging  from  the  aperture  is  aberration  free,  if  its 

wavefront  passing  through  the  center  of  the  aperture  is  spherical  with  a  center 

of  curvature  at  the  focus.  Any  deviations  of  the  wavefront  from  this  spherical 

form  represent  aberrations.  In  that  case,  the  spherical  wavefront  centered  at 

the  focus  is  referred  to  as  the  reference  sphere.  For  small  aberrations,  the 

Strehl  ratio  (i.e.,  the  ratio  of  central  irradiances  with  and  without  an 

aberration)  is  determined  by  the  aberration  variance  across  the 

24 

amplitude -weigh ted  aperture.  The  problem  of  aberration  balancing  for 

25  26 

uniformly  illuminated  circular  apertures  was  discussed  by  Nijboer.  *  This 

.  24  27 

has  been  extended  to  nonuniformly  illuminated  apertures.  *  In  aberration 
balancing,  a  classical  aberration  of  a  certain  order  (which  represents  a  term  in 
the  power  series  expansion  of  the  aberration  function  in  aperture  coordinates) 
is  mixed  with  aberrations  of  lower  order  such  that  the  variance  of  the  net 
aberration  is  minimized.  Consider,  for  example,  a  typical  balanced 
aberration2^  (representing  a  term  in  the  expansion  of  the  aberration  in  terms 
of  a  set  of  "Zernike"  polynomials  which  are  orthonormal  over  the  amplitude- 
weighted  annular  aperture) 

4>(P,0;e)  =  c  e  /2(n  +  1)  Rn’(p;e)cosm0,  (43) 

nm  m  n 

where  (p,0)  are  the  polar  coordinates  of  a  point  in  the  aperture  plane,  n 

and  m  are  positive  integers  (including  zero),  n  -  m  >  0  and  even,  Rm(p;e) 

—  n 

is  a  radial  polynomial  of  degree  n  in  p  and  has  the  form 


nm<  \  m  n  .  m  n— 2  m  m 

R  (p;e)  =ap+bp  +  ...  +  d  p 
n  n  n  n 


(44) 


2^*Mahajan,  V.  N. ,  "Zernike  Annular  Polynomials  for  Imaging  Systems  with 
Annular  Pupils,"  J.  Opt .  Soc .  Am. ,  71 ,  1981,  pp.  75-85,  1408. 

2^Nijboer,  B.  R.  A.,  'The  Diffraction  Theory  of  Aberrations,"  Ph.D.  Thesis 
(University  of  Groningen,  Groningen,  The  Netherlands,  1942). 

2^See  Ref.  18,  Chapter  9. 

27Szapiel,  S.,  "Aberration  Balancing  Technique  for  Radially  Symmetric 

Amplitude  Distributions:  A  Generalization  of  the  Marechal  Approach,"  J.  Opt . 
Soc .  Am. ,  72 ,  1982,  pp.  947-956. 
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The  coefficients  am  et  cetera,  depend  on  the  obscuration  ratio  t  and  the 
n 

aperture  amplitude  /iTpT.  The  quantity 


=  \ 

’/2  ’ 

m  =  0 

-  1, 

m  f  0 

(45) 

Unless  n  =  m  =  0,  the  coefficient  c  represents  the  standard  deviation  of 

nm 

the  aberration,  i.e. 


where 


and 


c  =  a. 

nm  9 


2  2  2 
=  <♦  >  ~  «&> 


1  2  it  /I  Ztt 

«t>">  =  f  J  A(p)  $n(p,e)pdpd ejf  J  /Up) 

£  0  ft  0 


pdpdB 


(46) 

(47) 
'  ( 

(48) 


From  Eqs.  (43)  and  (46),  we  note  that  the  standard  deviation  of  an  aberration 
can  be  obtained  immediately  by  comparing  its  form  with  the  corresponding 
orthonormal  aberration  represented  by  Eq.  (43),  without  having  to  calculate  the 
integrals  in  Eq.  (48).  The  variance  of  an  aberration  consisting  of  two  or  more 
terms  of  the  form  of  Eq.  (43)  is  given  by  the  sum  of  the  variance  of  each  of 
Che  aberration  terms. 


The  radial  polynomials  corresponding  to  balanced  primary  aberrations  are 
listed  in  Table  2  for  both  uniform  and  Gaussian  beams.  We  now  consider 
spherical  aberration,  coma,  and  astigmatism  in  uniform  and  Gaussian  beams  and 
compare  their  balancing  with  other  aberrations  to  minimize  their  variance 
across  the  aperture. 

6.1  SPHERICAL  ABERRATION 


First,  we  consider  spherical  aberration 

Vp)  =  v4 


(49) 


A 
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Table  2.  Radial  Polynomials  for  Balanced  Primary 

Aberrations  for  Uniform  and  Gaussian  Beams 


1 


-41- 


I 

f 


r 


where  A  represents  the  peak  value  of  the  aberration.  For  a  uniform  beam, 

S  .  .  ...  24 

the  standard  deviation  of  the  aberration  is  given  by 

a  =  [(4  -  e2  -  6e“  -  e6  +  4e8) 1/2/3/5]A  (50) 

su  s 

2 

If  we  mix  the  aberration  with  an  appropriate  amount  of  defocus  A^p  ,  we 
obtain  the  balanced  aberration 


*.  (  p)  =  A  p*  +  A.p2  (51  ) 

bs  s  d 

The  value  of  A,  which  minimizes  the  variance  of  the  balanced  aberration  can 
a 

be  obtained  by  comparing  it  with  the  radial  polynomial  R^(p;e)  given 
in  Table  2.  For  a  uniform  beam  it  is  given  by 


A  =  -  (  1  +  e2)A 
du  s 


(52) 


The  corresponding  standard  deviation  of  the  balanced  aberration  is  given  by2^ 


a  =  ((1  -  e2)  /6/?]A 
bsu  s 


(53) 


Note  that  for  a  given  value  of  A  ,  a,  decreases  as  e  increases, 

s  bsu 


For  a  Gaussian  circular  beam  with  y  =  1 

a  =  [ ( 20e 2  -  69e  +  40)1/2/(e  -  1)  ] A 
sg  s 

=  A  /  3 . 6  6  8 
s 


( 54a ) 
(54b) 


The  balancing  defocus  and  the  standard  deviation  of  the  corresponding  balanced 
aberration  are  given  by 


and 


dg 


(b°/a”)A 
4  4s 


-  0.933  A 


bsg 


=  A  //5  a? 
s  4 

=  A  / 1 3. 705 
s 


(55a) 

(55b) 

(56a) 

(56b) 
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Comparing  uniform  and  Gaussian  circular  beams  aberrated  by  spherical 
aberration,  we  note  the  following.  For  a  given  value  of 

a  /a  =  0.91  (5 

sg  su 

hence,  for  a  given  Strehl  ratio  exp(-a2),  where  a 2  is  the  variance  of  the 
28 

phase  aberration,  a  Gaussian  beam  can  tolerate  a  slightly  higher  amount  of 
spherical  aberration  than  a  uniform  beam.  By  balancing  with  an  appropriate 
amount  of  defocus,  the  standard  deviation  for  a  Gaussian  beam  is  reduced  by  a 
factor  of  3.74  compared  to  a  factor  of  4  for  a  uniform  beam.  Comparing  the 
standard  deviation  of  the  balanced  aberration  for  the  two  beams,  we  find  that 


O.  /a.  =0.98 

bsg  bsu 

Hence,  for  a  given  value  of  A  ,  the  Strehl  ratios  for  the  two  beams  are 

s 

practically  the  same.  We  noted  earlier  that  the  central  irradiance  for  an 
aberration-free  Gaussian  beam  is  lower  by  7.6%  compared  to  that  for  a 
corresponding  unifor-  beam.  For  beams  aberrated  with  a  small  value  of  A^ , 
the  difference  in  their  peak  central  irradiances  will  be  slightly  smaller. 
Comparing  the  amounts  of  balancing  defocus 


A.  /A  =  0.933  '  (59 

dg  du 

we  note  that  the  defocused  plane  for  a  Gaussian  beam  is  closer  to  the  focal 
plane  than  that  for  a  uniform  beam.  The  location  of  the  defocused  plane,  i.e., 
its  z  value,  is  given  by 

4>  ■  AJ  (60! 

o  d 


For  z  *  R,  it  is  given  by 


z  =  R  -  8AF  A 


where  A^  is  in  units  of  A. 
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To  compare  the  central  irradiance  in  the  focal  plane,  we  note  that  for  a 
quarter  wave  of  spherical  aberration,  a  Strehl  ratio  of  0.8  is  obtained  for  a 

uniform  beamio  correspond ing  to  au  =  0.22.  The  corresponding  Strehl  ratio 

-0. 22(0. 9 1 ) 2 

for  a  Gaussian  beam  is  e  ”  =  0.83,  giving  a  central  irradiance  of 

0.77  (compared  to  0.8  for  a  uniform  beam). 

6 . 2  COMA 


Next  we  consider  coma 


9  (p,0)  =  A  p3  cos© 
c  c 


(62) 


where  A  represents  the  peak  value  of  the  aberration.  For  a  uniform  beam, 

c  24 

the  standard  deviation  of  the  aberration  is  given  by 


a  *  [(1  +  e2  +  e"  +  e6)/8]l/2A 

cu  c 


(63) 


If  we  mix  coma  with  an  appropriate  amount  of  tilt,  Atp  cos8,  we  obtain 
the  balanced  aberration 


(p,0)  =  (A  p3  +  A  p)  cos0 
be  c  t 


(64) 


The  value  of  A  ,  which  minimizes  the  variance  of  the  balanced  aberration,  can 

_  | 

be  obtained  by  comparing  it  with  the  radial  polynomial  R^(p;e).  For  a 
uniform  beam,  it  is  given  by 


tu 


(2/3) [(1  +  t2  +  e“)/(l  +  e2) ]A 


The  corresponding  standard  deviation  of  the  balanced  aberration  is  given  by 


(65) 

24 


°b 


*  (1  ~  e2)(l  ♦  4e2  e1*)^2 

cu  *  6/2  (1  ^  e2)1/2 


(66) 


^^Mahajan,  V.  N. ,  "Strehl  Ratio  for  Primary  Aberrations  in  Terms  of  Their 
Aberration  Variance,"  J.  Opt .  Soc .  Am. ,  73 ,  1983,  pp.  860-861. 
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For  a  given  value  of  A  ,  <7.  slightly  increases  as  e  increases,  achieves  its 
®  c  bcu  °  ’ 

maximum  value  at  £  =  0.29,  and  decreases  rapidly  for  larger  values  of  e. 


For  a  Gaussian  circular  beam  with  y  =  1 


A  =  (b^/ahA 
tg  3  3c 


(67a) 


=  -0.608  A 


(67b) 


a,  =  A  /2/2  al 
beg  c  3 


=  A  /8.802 
c 


(68a) 


(68b) 


Comparing  uniform  and  Gaussian  circular  beams  aberrated  by  coma,  we  note 
the  following.  For  a  given  value  of  A^ 

■W°bc„  ’  °-96  <6! 

Therefore,  the  Strehl  ratio  for  a  Gaussian  beam  is  approximately  0.82  when  it 
is  0.80  for  a  uniform  beam.  Accordingly,  the  peak  irradiances  for  the  two 
beams  are  0.75  and  0.80,  respectively.  Note,  however,  that  since 

A  /A  =  0.91  ( 7( 

tg  tu 

the  peak  irradiance  for  a  Gaussian  beam  occurs  in  the  focal  plane  at  a  point 
which  is  closer  to  the  focal  point  by  9 X  compared  to  that  for  a  uniform  beam. 

6.3  ASTIGMATISM 

Finally,  we  consider  astigmatism 


$  ( p, 0)  =  A  o2  cos  8 
a  a 
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where  A  represents  the  peak  value  of  the  aberration.  For  a  uniform  beam, 
a 

the  standard  derivation  of  the  aberration  is  given  by 

a  =  ( 1/4) ( 1  +  e'*)1/2  A  (72) 

au  a 

If  we  mix  astigmatism  with  an  appropriate  amount  of  defocus,  we  obtain  the 
balanced  aberration 


4>  (p,0)  *  A  p2  cos20  +  A  p2  (73) 

ba  a  d 

For  both  uniform  and  Gaussian  beams,  the  amount  of  defocus  that  minimizes  the 

variance  of  the  balanced  aberration  can  be  obtained  by  comparing  it  with  the 

2 

Zernike  polynomial  R^( p; e)cos20.  It  is  given  by 

A  =  -  ( 1/ 2 ) A  (74) 

d  a 

so  that  the  balanced  aberration  becomes 

<t  ( p,  0)  =  ( 1/2) A  p2  cos20  (73) 

ba  a 

For  a  uniform  beam,  the  corresponding  standard  deviation  is  given  by 

a  =  ( 1/2/h) (1  +  e2  +  e-)1/2A  (76) 

bau  a 

For  a  given  value  of  A  ,  a  and  O,  both  increase  as  e  increases.  This 
®  a  au  bau 

is  true  for  a  Gaussian  beam  as  well.  For  a  Gaussian  circular  beam  with  y  -  1 


a  =  A  /2/6  a2  (77a) 

bag  a  2 


=  A  /5.609 
a 


(77b) 


Comparing  uniform  and  Gaussian  circular  beams,  we  note  that,  for  given  value 
of  A 

a 


°bag^°bau 


0.87 


(78) 
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Therefore,  the  Strehl  ratio  for  a  Gaussian  beam  is  approximately  0.85  compared 
to  a  value  of  0.80  for  a  uniform  beam.  Accordingly,  the  central  irradiances 
for  the  two  beams  are  0.78  and  0.80,  respectively. 

For  easy  comparison,  the  standard  deviation  of  primary  aberrations  and  of 
the  corresponding  balanced  aberrations  is  tabulated  in  Table  3  for  uniform  and 
Gaussian  (y  =  1)  circular  beams.  It  is  evident  that  the  standard  deviation 
for  a  Gaussian  beam  is  somewhat  smaller  than  the  corresponding  value  for  a 
uniform  beam.  Accordingly,  for  a  given  small  amount  of  aberration,  the  Strehl 
ratio  for  a  Gaussian  beam  is  higher  than  the  corresponding  value  for  a  uniform 
beam.  Similarly,  for  a  given  Strehl  ratio,  the  aberration  tolerance  for  a 
Gaussian  beam  is  somewhat  higher  than  that  for  a  uniform  beam.  We  also  note 
that,  whereas  aberration  balancing  in  the  case  of  a  uniform  beam  reduces  the 
standard  deviation  of  spherical  aberration  and  coma  by  factors  of  4  and  3, 
respectively,  the  reduction  in  the  case  of  astigmatism  is  only  a  factor  of 
1.22.  For  a  Gaussian  (y  =  1)  beam,  the  reduction  factors  have  a  similar 
trend,  but  the  factors  are  smaller:  3.74,  2.64,  and  1.16  corresponding  to 
spherical  aberration,  coma,  and  astigmatism,  respectively. 

We  note  that  the  variance  of  each  of  the  three  primary  aberrations  is 
minimized  if  it  is  measured  with  respect  to  a  reference  sphere  centered  at  a 
point  which  is  slightly  different  from  the  focal  point.  In  the  case  of 
spherical  aberration  and  astigmatism,  the  balancing  aberration  is  defocus, 
i.e.,  their  variance  is  minimized  if  they  are  measured  with  respect  to  a 
reference  sphere  centered  on  an  axial  point  z  f*  R.  In  the  case  of  coma, 
the  balancing  aberration  is  tilt;  i.e.,  its  variance  is  minimized  if  it  is 
measured  with  respect  to  a  reference  sphere  centered  at  a  point  in  the  focal 
plane  but  not  at  the  focal  point.  For  small  aberrations,  minimum  variance  has 
the  consequence  that  when  an  aberrated  beam  is  focused  on  a  target  at  a  fixed 
distance  z  from  the  aperture,  the  central  irradiance  on  it  is  maximum  if  the 
beam  is  slightly  defocused  in  the  case  of  spherical  aberration  and  astigmatism 
and  if  it  is  slightly  tilted  in  the  case  of  coma.  If  the  beam  is  neither 
defocused  nor  tilted,  then  we  can  say  that  for  focusing  systems  with  large 
Fresnel  numbers  (so  that  the  highest  aberration-free  irradiance  peak  occurs  at 
the  focal  point),  the  highest  peak  of  the  three-dimensional  irradiance 
distribution  occurs  at  the  point  with  respect  to  which  the  aberration  variance 
is  minimum. 
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Table  3.  Standard  Deviation  of  Primary  Aberrations  and  of  the  Corresponding 
Balanced  Aberrations  for  Uniform  and  Gaussian  (y  =  1) 

Circular  Beams 


Aberration,  $(p,9) 

Standard  Deviation 

Uniform  (y  =  0)  Gaussian  (y  =  l) 

Spherical 

Asp“ 

Ag/3. 35  As/3 . 67 

Balanced  Spherical 

Agp"*  +  Adp2 

As/13.42  Ag/13.71 

Coma 

Acp3cos0 

Ac/2.83  Ac/3.33 

Balanced  Coma 

(Acp3  +  Atp)cos6 

Ac/8.49  Ac/8.80 

Ast igmat ism 

Aap2cos28 

Aa/4  Aa/4.84 

Balanced  Astigmatism 

Aap2cos20  +  Adp2 

Aa/4.90  Aa/3.61 

“  (l/2)Aacos20 

- 

Adu  -  ~AS ,  Adg  *  -0.93A8  ;  Atu  =  — ( 2 / 3 ) Ac ,  Atg  =  -0.61AC  ; 

Adu  =  Adg  *  -Aa/2 
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7.  WEAKLY  TRUNCATED  GAUSSIAN  CIRCULAR  BEAMS 


7.1  IRRADIANCE  DISTRIBUTION  AND  BEAM  RADIUS 


In  this  section,  we  consider  a  weakly  truncated  Gaussian  circular  beam, 
i.e.,  one  for  which  the  aperture  radius  a  is  much  larger  than  the  beam  radius 
0).  When  a  »  oj,  i.e.,  when  y  is  very  large,  f(y;0)  ■+•  2y ;  therefore,  Eq.  (5) 
for  the  aperture  irradiance  distribution  may  be  written 

I(p)  =  2y  exp(-2yp2) 
or 

I(p)  =  (2P  /ttoj2)  exp[-2(p/oj)2] 
o 

In  Eq.  (79a),  I(p)  is  in  units  of  Pq/A  and  p  is  in  units  of  a  as  stipulated 
in  Eq.  (11).  In  Eq .  (79b),  no  such  normalization  is  used. 

For  large  y,  the  upper  limit  of  integration  in  Eq.  (11)  may  be  replaced 

by  infinity  with  a  negligible  error.  Since  we  are  considering  a  circular 

beam,  the  lower  limit  e  is  zero.  Hence,  if  we  let  B  =  irr  and  a  =  y  -  i*!1  , 

o 

the  integral  in  Eq.  (11)  reduces  to  the  form 


(79a) 

(79b) 


J  exp(-ap2) Jg(Bp)pdp  (80a) 


29 

which  for  Rea  >  0  is  equal  to 

(l/2a)  exp(-62/4a)  (80b) 

Accordingly,  it  can  be  shown  that,  for  large  y,  Eq.  ill)  reduces  to 

l(r;z)  =  (R/z)2f2y/(4>2  +  y2)  ]exp[-yir2r  2/2(<f2  +  y2)]  (81a) 


29cradshteyn,  I.  S.,  and  I.  M.  Ryzhik,  Tables  of  Integrals,  Series,  and 
Products ,  4th  ed.,  Academic,  NY,  1965,  p.  717. 


or 


(81b) 


I(r;z)  =  (2P  /mo2)  exp(-2r 2/<o2) 
0  2  2 


where 


to2  =  (Xz/mo)2  +  to2(l  -  z/R)2  (82) 

z 

In  Eq.  (81a),  the  irradiance  is  in  units  of  P  A/X2R2  and  r  is  in  units  of 

o 

Xz/D  as  was  the  case  in  Eq.  (11).  In  Eq.  (81b),  these  quantities  are  not 
normalized . 

Comparing  Eqs.  (79)  and  (81),  we  note  the  well-known  result  that,  if  the 
truncation  of  the  beam  by  the  aperture  is  neglected,  a  Gaussian  beam  propagates 
as  a  Gaussian  beam.  The  beam  radius  at  a  distance  z  from  the  aperture  is  to 

2 

given  by  Eq .  (82). 

If  we  let  r  -  0  in  Eq.  (81),  we  obtain  the  axial  irradiance 

I ( 0 ; z )  =  (R/z;2[2y/($2  +  y2)]  (83a) 

o 

or 

I ( 0 ; z )  =  2P  /mo2  (83b) 

o  z 

Of  course,  for  large  y  and  e  =  0,  Eq.  (24)  reduces  to  Eq.  (83)  as  expected. 

If  we  let  z  =  R  in  Eq.  (81),  we  obtain  the  focal-plane  irradiance 
distribution 


I(r;R)  =  (2/y)exp(-TT2r2/2y) 


(84a) 


or 

I(r ;R)  *  (2P  /mu2)  exp(-2r2/co2 ) 

OK  K 


where 


to  *  XR/mu 
R 


(84b) 


(84c) 
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is  the  beam  radius  in  the  focal  plane.  The  focal-point  irradiance  is  given  by 


1(0 ;R)  =  2/y  (85a) 

I(0;R)  =  2ttP  u)2/\2R2  (85b) 

o 

a  result  already  obtained  in  Eq.  (22). 

If  we  equate  to  zero  the  derivative  of  Eq.  (83)  with  respect  to  z,  we 
obtain  the  z  value  at  which  the  axial  irradiance  is  maximum.  Note  that  this 
value  is  given  by 


z  /R  =  [1  +  (y/irN)2]-1 
P 


(86a) 


or 


z  /R  =  (1  +  (AR/irur)  ] 
P 


(86b) 


Substituting  Eq.  (86)  into  Eq.  (83),  we  obtair  the  peak  value  of  axial  irradiance 


I ( 0 ; z  )  =  (2/y)  +  (2y/irzN2) 
P 


(87a) 


l(0;z  )  =  2P  / m j 

p  o  zp 


(87b) 


where 


oj2  =  u)2/[l  +  (ttu)2/AR)2] 

zp 


(87c) 


Comparing  Eqs.  (85a)  and  (87a),  we  note  that  the  peak  axial  irradiance  is  higher 

than  the  corresponding  focal-point  irradiance  by  2y/iT2N2  (in  units  of  P  A/A2R2) 

o 

or  by  2Po/ttw2.  Equations  (86a)  and  (87c)  can  also  be  written  in  the  form 


z  /R  =  [1  +  (ttN  )~2]-1 
P  g 


(86c) 


and 


oj2  -  a)2/ [  1  +  (ttN  )2] 
zp  g 


(87d) 


-51- 


i  12 

respectively,  where  N  =  u>  /XR  is  the  Gaussian  Fresnel  number.  It  repre- 

§ 

sents  the  number  of  Fresnel  zones  in  the  aperture  plane  within  the  Gaussian 

beam  radius  as  observed  from  the  focus,  just  as  N  represents  the  number  of 

zones  within  the  full  aperture,  provided  that  there  is  no  obscuration.  Since 

the  axial  irradiance  is  maximum  at  z  ,  the  beam  radius  is  minimum  at  this 

P 

position.  This  may  also  be  seen  by  equating  to  zero  the  derivative  of  Eq.  (82) 

with  respect  to  z.  The  minimum  beam  radius  u)  is  generally  referred  to  in  the 

.1  ZP 

literature  as  the  beam  waist. 

It  is  evident  that  z  <  R,  i.e.,  the  peak  of  axial  irradiance  does  not 

P 

occur  at  the  focal  point  but  at  a  point  between  it  and  the  aperture.  Moreover, 

u)  <  a),  i.e.,  the  waist  of  the  diffracted  beam  is  smaller  than  the  beam  radius 
zp 

at  the  aperture.  Note,  however,  that,  as  discussed  in  Section  4,  even  though 
the  peak  axial  irradiance  and  waist  are  not  located  at  the  focal  point  (z  =  R) , 
the  smallest  beam  radius  and  maximum  central  irradiance  on  a  target  at  a  fixed 
distance  z  are  obtained  when  the  beam  is  focused  on  it.* 


The  encircled-power  distribution  in  an  observation  plane  is  given  by 

P(r  ;z)  =  1  -  exp[-Yir2r2/2(<I>2  +  y^)  ]  (88a) 

o  o  o 

or 

P(r  ;z)  =  P  (1  -  exp[-2r2/u)2)  ]  (88b) 

o  o  o  z 

where  r  is  in  units  of  Xz/D  in  Eq .  (88a). 
o  ^ 


The  radius  of 
z  from  the  aperture 
in  Eq.  (1).  Noting 


curvature  R  of  the  Gaussian  spherical  wave  at  a  distance 
z 

is  given  by  the  inverse  of  the  coefficient  of  -iur2/X 
that  r  is  in  units  of  Xz/D  in  8  =  tt r ,  we  find  that 


z/R  *  (A/Xz)lm(l/a)  -  1 
z 

-  I A4>  /Xz<4>2  ♦  y2)]  -  1 
o  o 


(89a) 


_ 1  -  z/R _ 

(1  -  z/R)2  +  ( X z / TTU) 2 )  2 


(89b) 


4 
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At  the  waist 
focal  plane, 
wave . 


position  z  ,  R  =  °°,  implying  a  plane  wave.  Similarly,  at  the 
p  zp 

Rz  =  -R.  A  negative  value  of  Rz  indicates  a  diverging  spherical 


The  results  for  a  weakly  truncated  collimated  Gaussian  beam  can  be 
obtained  from  those  for  a  focused  beam  by  letting  R  -*■  00 .  Thus,  for  example 

I(r;z)  =  {2y/[l  +  (4yz/tt)2]}  exp{-8Yz2r2/[l  +  (4yz/tt)2]}  (90a) 

or 

I(r;z)  =  (2P  /™2)  exp(-2r 2/iij2)  (90b) 

o  z  z 

where 

U)2  =  to2[l  +  (Xz/Tto2)2]  (90c) 

z 

and 

R  =  -z[l  +  (  TToj  2  /  X  z  )  2  ]  (90d  ) 

z 

In  Eq.  (90a),  the  irradiance  is  in  units  of  P  /A,  r  is  in  units  of  Az/D 

^  2 

(z  is  not  normalized  here),  and  z  is  in  units  of  D  /X.  We  note  from  Eq .  (90c) 

that  to  increases  monotonically  as  z  increases,  i.e.,  the  beam  expands  as  it 
z 

propagates . 

In  this  section,  we  have  written  equations  in  two  equivalent  forms. 
Equations  (a)  are  written  in  a  normalized  form  so  that  they  can  be  investigated 
parametrically.  Equations  (b)  are  written  without  any  normal izat ions ,  and  they 
are  more  suitable  for  evaluating  results  when  the  specific  parameters  involved 
are  known.  Incidentally,  note  that  all  equations  in  this  section  reduce  to  the 
corresponding  quantities  in  the  aperture  plane  if  we  let  z  0. 

Figure  9  shows  the  focal-plane  irradiance  and  encircled-power  distribu¬ 
tions  for  /y  *  2.  The  solid  curves  have  been  obtained  by  using  Eqs.  (17)  and 
(18)  [or  Eqs.  (30  and  (32)  with  <t>o  =0],  and  the  dashed  curves  represent 
their  corresponding  approximations  given  by  Eqs.  (81a)  and  (88a),  respectively, 
with  *  0.  We  note  that  the  approximate  results  agree  well  with  the  true 
results.  The  maximum  difference,  which  occurs  at  the  focus,  is  less  than  4%. 

For  larger  Y»  the  agreement  is  found  to  be  even  better. 
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Note:  The  irradiance  and  encircled  power  are  in  units  of  P0A/A2R2  and 

P0,  respectively.  The  radial  distance  r  or  rG  in  the  focal  plane  are 
in  units  of  AR/D.  The  focal  point  is  at  r  "  0. 


Figure  9.  Focal-Plane  Irradiance  and  Enc ire led-Power  Distributions  for  a 
Gaussian  Beam  with  /y  ■  2 


I 

4 


Figure  10  shows  how  the  axial  irradiance  of  a  focused  Gaussian  beam  varies 

when  /y  =  2  and  N  =  1,  10,  and  100.  The  irradiance  is  in  units  of  P  A/A2R2. 

o 

Once  again,  the  solid  curves  in  this  figure  have  been  obtained  by  using  Eq  . 

(24),  and  the  dashed  curves  represent  their  corresponding  approximations  given 
by  Eq.  (83a).  It  is  evident  that  Eq.  (83a)  represents  the  true  axial  irradiance 
quite  well.  The  only  significant  difference  occurs  when  N  =  1,  in  that  the  true 
results  show  secondary  maxima  and  minima  but  the  approximate  result  shows  only 
the  principal  maximum.  For  larger  values  of  y,  e.g.,  /y  =  2.5,  the 
secondary  maxima  and  minima  disappear,  and  the  true  and  approximate  results 
overlap  each  other  at  the  scale  of  Figure  10.  Hence,  we  conclude  that  the 
truncation  of  an  aberration-free  Gaussian  beam  by  an  aperture  has  a  negligible 
effect  on  the  irradiance  distribution  as  the  beam  propagates  when  /y  >_  2. 

7.2  ABERRATION  BALANCING 


When  a  Gaussian  circular  beam  is  weakly  truncated,  i.e.,  when  y  is 
large,  the  quantity  pg  in  Table  2  reduces  to 

ps  =  <ps> 

=  (s/2y)  pg_2 

=  (s/2) :  y_s/ 2  (91) 

As  a  result,  we  obtain  simple  expressions  for  the  radial  polynomials  represent¬ 
ing  balanced  primary  aberrations.  They  are  listed  in  Table  2.  If  we  normalize 
0  by  u  (instead  of  by  a  as  in  Table  2),  then  y  disappears  from  these 
expressions.  As  in  Section  6,  the  standard  deviation  of  an  aberration  can  be 
obtained  by  comparing  its  iOrm  with  the  corresponding  orthonormal  aberration  of 
Eq.  (43). 

The  standard  deviation  of  primary  aberrations  and  of  the  corresponding 
balanced  aberrations  for  a  weakly  truncated  Gaussian  beam  is  given  in  Table  4. 

In  this  table,  p'  a  /y  p  is  a  radial  variable  in  the  aperture  plane  normalized 
by  the  beam  radius  U),  and  the  aberration  coefficients  A|  represent  the  peak 
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100 


value  of  a  classical  aberration  at  p'  =  1.  The  reason  for  defining  the 
aberration  coefficients  in  this  manner  should  be  obvious.  Since  the  power  in  a 
very  weakly  truncated  Gaussian  beam  is  concentrated  in  a  small  region  near  the 
center  of  the  aperture,  the  effect  of  the  aberration  in  its  outer  region  is 
negligible.  Accordingly,  aberration  tolerances  in  terms  of  the  peak  value  of 
the  aberration  at  the  edge  of  the  aperture  are  not  very  meaningful. 

We  note  from  Table  2  (or  Table  4)  that  the  point  with  respect  to  which  the 
variance  of  an  aberration  is  minimized  is  given  by 

$  =  ~(4/y)A  (92a) 

o  s 

=  -4y  A’ 

1  s 

A  =  -( 2/y) A  (92b) 

t  c 

2/y  A' 
c 

( 1/2)A 

a 

(  y/2) A^  (92c) 

in  the  case  of  spherical  aberration,  coma,  and  astigmatism,  respectively.  From 

Table  4,  we  note  that  the  balancing  of  a  primary  aberration  reduces  its  standard 

deviation  by  a  factor  of  /5 ,  /3 ,  and  /2  in  the  case  of  spherical  aberration, 

coma,  and  astigmatism,  respectively.  These  reduction  factors  are  listed  in 

_2 

Table  5  for  uniform  (y  =  0)  beam,  and  e  truncated  (y  =  1)  and  weakly 
truncated  (large  y)  Gaussian  beams.  It  is  evident  that  as  y  increases,  the 
reduction  factors  decrease  for  each  of  the  three  primary  aberrations.  The 
amount  of  balancing  aberration  decreases  as  y  increases  in  the  case  of 
spherical  aberration  and  coma,  but  it  does  not  change  in  the  case  of  astig¬ 
matism.  For  example,  in  the  case  of  spherical  aberration,  the  amount  of 
balancing  defocus  for  a  weakly  truncated  Gaussian  beam  is  (4/y)  times  the 
corresponding  amount  for  a  uniform  beam.  Similarly,  in  the  case  of  coma,  the 
balancing  tilt  for  a  weakly  truncated  Gaussian  beam  is  (3/y)  times  the 
corresponding  amount  for  a  uniform  beam. 


and 
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Table  4.  Standard  Deviation  of  Primary  Aberrations  in  Weakly 
Truncated  (/y" >  3)  Gaussian  Circular  Beams 


Aberration,  4>(p,0) 

a* 

Spherical 

A'p’  “ 

2/5  A' 

s 

s 

Balanced  Spherical 

A’(p'4  -  4p'2) 

2  A' 

s 

s 

Coma 

A'  p'  3  cos  0 

/3  A' 

c 

c 

Balanced  Coma 

A' ( p' 3  -  2 p' )cos 0 

A' 

c 

c 

Astigmatism 

A'p'  2cos20 
a 

(1//2)  A' 
a 

Balanced  Astigmatism 

A'p' 2cos20  -  ( 1/2)A' p' 2 

(1/2)  A' 

a  a 

a 

=  (1/2)  A’  p'  2  cos2 0 

a 

p'  =  /yp,  Ag  -  Ag/y2,  A£  *  Ac/y3/2 

,  Ag  =  Aa/y,  A^  =  Ad/y,  A£  =  At//y. 

Table  5.  Standard  Deviation  Reduction  Factor  for  Primary  Aberrations* 


Ah er rat i on 

Reduction  Factor 

Uniform  ( Y  =  0) 

Gaussian  ( y  *  1 ) 

Weaklv  Truncated 
Gaussian  (»/y  >  3) 

Spherical 

4 

3.74 

2.24 

Coma 

3 

2.A4 

1.73 

Astigmatism 

1.22 

1.A6 

1 .41 

*It  represents  the  factor  bv  which  the  standard  deviation  of  a  classical 
aberration  across  a  circular  aperture  is  reduced  when  it  is  optimally  balanced 
with  other  aberrations. 


Table  6  gives  the  reduction  factors  that  relate  the  peak  value  A^  of  a 
primary  aberration  at  the  edge  of  a  circular  aperture  and  the  standard  devia¬ 
tion  of  its  corresponding  balanced  aberration  for  various  values  of  y.  In  the 
case  of  balanced  aberrations,  these  numbers  are  given  by  /F  a^,  2/2  a^,  and 
2/6  a^  for  spherical  aberration,  coma,  and  astigmatism,  respectively.  For  exam¬ 
ple,  for  spherical  aberration  A  p4,  the  standard  deviation  of  the  corresponding 

s 

balanced  spherical  aberration  when  /y  '  2  is  equal  to  A^/18.29. 

Comparing  the  standard  deviation  results  given  in  Table  6  with  those  for  a 

weakly  truncated  Gaussian  beam  in  Table  4,  we  find  that  they  agree  with  each 

other  with  negligible  difference  for  /y  >  3.  This  provides  a  convenient 

definition  for  a  weakly  truncated  Gaussian  beam,  namely,  that  a  >  3u). 

30  3 1 

Some  authors  ’  have  assumed  that  /y  ^  2  provides  a  sufficient  condition 

for  the  validity  of  the  aberration  analysis  of  a  weakly  truncated  Gaussian 

beam  given  here.  When  /y  =  2 ,  the  standard  deviation  of  balanced  spherical 

aberration  according  to  the  weakly  truncated  beam  assumption  is  given  by 

A  /8,  while  the  true  value,  as  stated  above,  is  given  by  A  /18.29,  which  is 
s  s 

significantly  different.  When  /y  =  3,  the  corresponding  standard  deviations 

are  given  by  A  /40.50  and  A  /43.52,  which  are  nearly  equal  to  each  other, 
s  s 

The  difference  between  the  true  and  approximate  results  is  even  less  for 
/y  >  3. 


When  /y  =  2,  even  though  the  true  focal-plane  distribution  obtained  from 
Eq.  (17)  or  Eq .  (30)  agrees  quite  well  with  that  obtained  from  Eq.  (84a),  the 
true  and  approximate  standard  deviations  of  primary  aberrations  are 
significantly  different  as  pointed  out  above.  The  reason  for  the  discrepancy  in 
the  case  of  an  aberrated  beam  is  simple.  Even  though  the  irradiance  in  the 
region  of  the  aperture  to/a  ^  p  <_  1  is  quite  small  compared  to  that  at  or 
near  its  center,  the  amplitude  in  this  region  is  not  as  small.  Moreover,  the 
aberration  in  this  region  can  be  quite  large  and  thus  have  a  significant  effect 
on  the  standard  deviation.  In  the  case  of  spherical  aberration,  it  increases 


^Yoshida,  A.,  "Spherical  Aberration  in  Ream  Optical  Systems,"  Appl.  Opt.,  21 
1982,  pp.  1812-1816.  -  “ 

■^Herloski.  R.,  "Strehl  Ratio  for  Untruncated  Aberrated  Gaussian  Beams," 

J.  Opt,  Soc.  Am..  A2,  1985,  pp.  1027-1030. 
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Table  6.  Standard  Deviation  Factor  for  Primary  Aberrations  for  a  Gaussian 
Circular  Beam  with  Various  Values  of  y 


/y 

Balanced 

Spherical 

Balanced 

Coma 

Balanced 

Astigmatism 

0 

13.42 

8.49 

4.90 

0.5 

13.69 

8.53 

5.06 

1.0 

13.71 

8.80 

5.61 

1.5 

14.90 

9.74 

6.81 

2.0 

18.29 

12.21 

9.08 

2.5 

26.33 

17.62 

12.82 

3.0 

43.52 

27.57 

18.06 

3.5 

75.78 

42.96 

24.51 

4.0 

128.09 

64.01 

32.00 

Note : 


The  numbers  represent  the  factor  by  which  the  peak  aberration  coefficient 
must  be  divided  by  in  order  to  obtain  the  standard  deviation. 
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as  p1*.  In  the  case  of  coma  and  astigmatism,  it  increases  as  p3  and  p2, 
respectively.  Hence,  we  require  a  larger  value  of  y,  namely,  /y  >  3, 
for  the  aberrated-beam  analysis  of  this  section  to  be  valid.  This  is  also 
true  of  defocus  which  varies  as  p2. 
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8.  CONCLUSIONS 


We  have  compared  the  effects  of  diffraction,  obscuration,  and  aberrations 

32 

on  the  propagation  of  uniform  and  Gaussian  beams.  This  comparison  is  based 
on  fixed  total  power  in  the  beam  regardless  of  the  obscuration  or  the  amplitude 
across  it  at  the  aperture  plane.  The  following  general  conclusions  can  be  drawn 
from  the  work  reported  here. 

The  mirrors  in  a  high-power  Gaussian  beam  train  are  illuminated  unevenly 
and  must  withstand  higher  flux  densities  than  those  in  a  corresponding  uniform 
beam  train.  The  uneven  illumination  will  cause  thermal  distortions  of  the 
mirrors  thereby  introducing  aberrations  in  the  beam. 

The  focal-point  irradiance  for  a  focused  Gaussian  beam  is  smaller  than  the 
corresponding  value  for  a  uniform  beam  of  the  same  total  power.  Also,  the 
encircled  power  for  small  circles  is  higher  for  a  uniform  circular  beam,  but  for 
large  circles,  it  is  higher  for  a  Gaussian  circular  beam. 

The  effect  of  the  Gaussian  apodization  of  the  aperture  is  to  Increase  the 
size  of  the  central  disc  and  to  decrease  the  power  in  the  rings  of  the 
diffraction  pattern.  The  effect  of  a  central  obscuration  Is  opposite  to  that  of 
the  Gaussian  apodization.  It  reduces  the  size  of  the  central  disc  and  increases 
the  power  in  the  rings  of  the  diffraction  pattern.  Accordingly,  as  the 
obscuration  Increases,  the  difference  between  the  diffraction  effects  of  uniform 
and  Gaussian  beams  decreases. 

The  minima  of  axial  irradiance  for  a  uniform  beam  have  a  value  of  zero, 
while  those  for  a  Gaussian  beam  have  nonzero  values.  While  the  principal 
maximum  of  axial  irradiance  for  a  Gaussian  beam  has  a  smaller  value  than  the 
corresponding  value  for  a  uniform  beam,  the  secondary  maxima  for  a  Gaussian  beam 
have  higher  values.  Even  though  the  principal  maximum  does  not  necessarily 


^Mahajan,  V.  N. ,  "Comparison  of  Uniform  and  Gaussian  Beam  Diffraction," 
SPIE  Proc.,  Diffraction  Phenomena  in  Optical  Engineering  Applications, 
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it,  unlike  the  conclusions  of  some  authors. 

For  a  <  u),  the  Gaussian  beams  are  somewhat  less  sensitive  to  aberrations 
than  uniform  beams.  Accordingly,  aberration  tolerance  is  somewhat  higher  for 
the  Gaussian  beams.  However,  this  tolerance  increases  rapidly  as  a  becomes 
larger  and  larger  compared  to  ui.  This  is  understandable  since,  for  a  much 
greater  than  u),  the  power  in  the  aperture  is  concentrated  in  a  small  region 
near  the  center  of  the  aperture;  therefore,  the  aberration  in  its  outer  region 
has  little  effect  on  the  irradiance  distribution.  For  a  >_  2io,  the  truncation 
of  an  aberration-free  Gaussian  beam  by  the  aperture  has  a  negligible  effect  on 
its  propagation.  Accordingly,  it  remains  a  Gaussian  beam  as  it  propagates. 
However,  when  the  beam  is  aberrated,  a  3<i)  is  required  in  order  to  neglect  the 
effect  of  its  truncation  by  the  aperture.  When  a  >_  3d),  it  is  more  appropriate 
to  define  aberration  coefficients  as  the  peak  aberrations  at  the  beam  radius 
ui  rather  than  at  the  aperture  edge  a,  since  the  power  in  the  beam  is 
concentrated  in  a  small  region  near  the  center  of  the  aperture,  and  the  effect 
of  an  aberration  in  its  outer  region  is  negligible.  With  the  aberration 
coefficients  defined  in  this  manner,  the  beam  becomes  most  sensitive  to 
spherical  aberration  and  least  sensitive  to  astigmatism,  rather  than  being  most 
sensitive  to  coma  and  least  sensitive  to  spherical  aberration. 


^Tanaka,  K. ,  N.  Saga,  and  K.  Hauchi,  "Focusing  of  a  Gaussian  Beam  Through  a 
Finite  Aperture  Lens,"  Appl .  Opt . ,  24,  1985,  pp.  1098-1101. 
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